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ABSTRACT 

We present a new technique to quantify the Hght contribution coming from the faint high 
redshift (z ^ 6) galaxies below the detection threshold of imaging data, set conventionally 
at S/N — 4.5. We illustrate the technique with an application to Hubble Space Telescope 
Advanced Camera for Surveys images in the F775W and F850LP filters of the Ultra Deep 
Field parallel field NICP12. The aim of this analysis is to extend by a few magnitudes the faint 
end of the luminosity function at 2; ~ 6. After masking all the detected sources in the field 
we apply a Fast Fourier Transform to obtain the spatial power spectrum of the background 
signal. The power spectrum permits us to separate the background noise signal, the residuals 
due to the data reduction of the wide field, and the overall signal produced by faint galaxies. 
The ratio of the signal in the ^775 and 2:350 bands is used to estimate the contribution of the 
faint i-dropout objects. We rely on extensive Monte Carlo simulations to characterize various 
sources of uncertainty and quantify the number of faint i-dropout galaxies in the field. The 
analysis allows us to put constraints on the luminosity function at z '^ 6 down to Z850= 30 
mag, 2.5 mag fainter than with standard techniques on the same data. The data are consistent 
with a faint end slope of the luminosity function of a = —1.9. Assuming a specific set of 
values for the clumping factor, escape fraction, and spectral energy distribution, we find that 
the z ^ Q undetected galaxies down to zgso — 30 mag could have driven cosmic reionization. 

Key words: cosmology: early universe - galaxies: high-redshift 



1 INTRODUCTION 

In the last few years we have begun to investigate the cosmic reion- 
ization of hydrogen by using a variety of techniq ues. Reionization 
was c ompleted between 2 ~ 10 and z ~ 6 JLoeb & Barkanal 
l200lh . bracketed by the detection of Gui m-Peterson troughs in the 
spectra of z ~ 6 quasars JFan et al.l2006r) and by the Compton opti- 
cal depth measured by the Wilkinson Microwave Anisotropy Probe 
(WMAP; Jarosiketal. 2011). However, the nature of the sources 
responsible for cosmic reionization is still an open question. 

The UV photons expected to contribute to the reionization 
process come from mainly two sources, star forming galaxies and 
quasars. The former probably were fundamental for the hydro- 
gen reioniza tion, whereas the latter l ikely dominated th e helium 
reionization JShapiro & Giroux|[T987l iMadau et al.lll999h . Accret- 
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ing black holes in the redshift range between 2: = 6 and 2 = 8 
are rare and possibly obscured by significant amounts of gas and 
dust, so their UV emission probabl y does not contribu te signif- 
icantly toJhereionization_grocessJTreistereLa^ 201ln. Accord- 
ing to iFan et alj ( 1200 ll) , iMeiksid J2005l) , and IShankar & MathuJ 
( 120070 ■ 2 ~ 6 quasars are not able to provide enough ionizing 
photons to contribute significantly to reionization. 

Population III stars are also unlikely to dominate the ionizing 
UV production. In the standard scenario, in which only a single 
metal free star per dark halo is formed, the contribution to reioniza- 
tion given by Population III stars is only indirect. In fact, while the 
cumulative number of ionizing photons produced by them is not 
enough, they enrich the intergalactic medium of metals and allow 
the mor e efficient formation of fo llowing generations of Population 
II stars dTrenti & Stiavelli 2009). The se stars are likely located in 
faint galaxies dLoeb & Barkanall200lh . 



A large optical depth to electron scattering, consistent with 
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the Universe b eing reionized at z ^ 15, was found by WMAP 
JBennett et^ |2003 ; Komatsu et alj 1201 ll) . Since there are still 
strong u ncertainties on several key p arame ters concerning reion- 
ization, iKuhlen & Faucher-Giguerd 1 120121) stated that it is not 
clear if galaxies could have completed the reionization by z '^ 
6 and account for the large optical depth to electron scattering 
found by WMAP-7. On the other hand, most of the studies car- 
ried out in the last decade are focused on the critical star for- 
mation rate density (SFRD ). The Lyman break g alaxies detected 
at z ~ 6 jst iavelli et aT] l2004bl) and z y 7 JBouwens et alJ 
|2004 lBolton& Haehneh toOTi ; lOesch et al.l l2009) in the Hubble 
Ultra Deep Field (hereafter HUDF) are not bright enough to ac- 
count for the required ionizing budget. However, more recently 
iFinkelstein et alj ( 120 12h have found that the detected galaxies at 
2 ~ 6 are s ufficient to fu l ly re ionize the IGM if C'/f esc ~ 10. 
According to lBunker et al.l ( l2010l) and lFinkelstein et al.l ( l2010h . the 
observed z ^ 7 galaxies could have completed the reionization 
process if fesc ^ 50%. The observed faint-end slopes imply large 
numbers of undetected faint galaxies which likely provide the bulk 
of ioniz ing pho tons and played an important r ole in cosmic reion- 
ization JYan & Windhorst 2004; Richard et al. 20081: JMcLure et al.] 
bOld l lTrentTet al.„201Q ; ,Bouwens et al.,,2011,) . 

Unfortunately, direct probing of the faint galaxy population 
is outside present capabilities and will need to wait for the James 
Webb Space Telescope (JWST). Thus, in order to make further 
progress in this field before the launch of JWST, one needs to find 
ways to study the faint end of the galaxy luminosity function at 
z ^ 6 using the integrated light contribution of individually unde- 
tected galaxies. This idea is not new and it was tried in the context 
of high-z studies using the Spitzer Infrared Arra y Camera (IRAC) 
and A dvanced Camera for Surveys (ACS) data dKashlinskv et al.l 
uOOTT) and is, in fact, also analogous to the study of surface bright- 
ness fluctuations of nearby galaxies (Tonry & Schneider 1988). 

The aim of this paper is to describe a new technique to an- 
alyze the contribution of faint, individually undetected, galaxies 
where we indicate as undetected galaxies those galaxies charac- 
terized b^_beingjmnterthan the typical signal-to-noise (S/N) ratio 
of 4.5 dOesch et al.ll2009l) . We focus mostly on galaxies just below 
the detection threshold as such galaxies are likely to exist, given 
that the luminosity function (LF) of galaxies should not know the 
HU DF magnitude limit, and they might be sufficient for reioniza- 
tion JYan & Windhorstll2004l) . 

These galaxies are too faint to be measured individually, but 
their signal can be quantified using a Fourier Transform based tech- 
nique. The main difference here over the previous attempts at high 
redshift is our special attention to systematic effects. 

The Lyman break technique permits to identify star-forming 
galaxies at high redshift using color criteria derived from a multi- 
band photo metry in the region acros s the Lyman-cont inuum dis- 
continuity JSteidel & HamiltonI [19931 ; iGiavaliscdl 19981) . Compar- 
ing the flux of a galaxy within a set of two or three ad-hoc filters 
depending on the redshift, it can be noticed that no flux is collected 
by the bluer filter due to the Lyman a absorption, while the object 
is detected in the redder bands. 

In this paper we focus on i-dropout galaxies at z ~ 6 that 
can be studied using images obtained with the same camera and 
characterized by, essentially, the same systematic effects. Galaxies 
at 2 ~ 6 are expected to be undetected in the ACS F775W band 
(hereafter 2775) image and detected in the ACS F850LP (hereafter 
zgso) one, so they are named 1775 -dropout galaxies. 

The faint galaxy population we seek is detected by an excess 



Table 1. Properties of the NICP12 data. 



Filter 



HST 
orbits" 



Exposure 
Time" [s] 



Zero Point 
[AB mag] 



10 o- Magnitude" 
[AB mag] 



«775 
-2850 



23 

69 



54,000 
168,000 



25.654 
24.862 



28.44 
28.47 



" From lOesch et all J2007l) 
*> From lde Marchi et all booi) 

NOTES: Col. 1: filter name. Col. 2: number of HST orbits. Col. 3: to- 
tal exposure time in seconds. Col. 4: zero point in AB magnitude. Col 5: 
10 a limiting magnitudes in AB magnitude computed within apertures of 
O.lS'.'radius. 



in the power spectrum in the zgso image with respect to that of the 
2775 one, i.e., it is a Lyman break population. 

In Section 2, we describe the data used. In Section 3, we ex- 
plain the data analysis and we present the main steps of the power 
spectrum technique adopted to obtain the signal coming from the 
undetected galaxies, which are bona fide 2; ~ 6 candidates. In 
Section 4, we describe how to estimate the contribution to cosmic 
reionization using Monte Carlo simulations, we discuss our results 
and their implications for the reionization epoch. 



2 DATA SET 

To date, the deepest optical data available are from the HUDF 
project ( JBeckwith et al.l|20 06) and the Hubble Ultra Deep Field 
Follow-Up (hereafter UDF05; Oesch et al. 2007). These ultra deep 
observations of multiple fields have been obtained with ACS. In 
this first application of our technique we selected the parallel field 
NICP12 rather than the main HUDF field. The reason for this is 
that the NICP12 images were produced with a more advanced dat a 
reduction than the one used in the main field ( lOesch et al .1120070 . 
In particular, an herringbone effect introduced inadvertently by the 
subtraction of a combined, high S/N , bias frame (superbias) made 
from compressed bias frames containing a small electronic noise 
component was eliminated as well as electronic ghost images of 
the bright sources. Thus, the improved version of the NICP12 im- 
ages (v2.0) is the cleanest available. We are currently reprocessing 
the main HUDF field so as to bring it to the same standard and 
when this field will be available we plan to apply this analysis to it 
as well. 

The coordinates of the center of the NICP12 field are 



R.A.=03''33™03.60%Dec.=-27°41'01. 



(J2000.0). The char- 



acteristics of the observations for the 1775 and zgso images are spec- 
ified in Table|2] Throughout this paper we used the AB mag magni- 
tude system. The original images are characterized by a pixel scale 
of 30 mas pixeP^ and a side of 10500 pixels. They were rotated 
to have North up and East left. 



3 DATA ANALYSIS 

A first step to reveal the contribution of undetected galaxies was to 
mask all the light coming from stars and detected galaxies, includ- 
ing their bright halos. Once the appropriate mask had been created 
we applied the power spectrum technique to analyze the different 
contributions existing in the background. 
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3.1 Creation of the Mask 

Since the MULTIDRIZZLE task fKoeke moer et al J 20021) works on 
dithered observations of the same field and only a few pointings in 
the outer regions of the field can be combined, we used the weight 
maps associated to the v2.0 data to select the central region for the 
following analysis. The data we focused on were binned to have 
a scale of 90 mas pixeP^ and a size of 3500 x 3500 pixels. This 
rebinning is intended to limit the computational volume, but it does 
not affect the final result, as discussed in Section [3^ 

First of all, we had to remove all the defects and objects that 
could interfere with our analysis, such as the detected galaxies, 
their bright halos, residual cosmic rays, and bad pixels. Galaxies 
were detected independently in the {775 and zgso images usin g the 
SExtractor photometry package by ISertin & AmoutsI (119961) , ver- 
sion 2.5.0. The SExtractor parameters were optimized to maximize 
the number of detected galaxies while minimizing the number of 
spurious sources. The detection threshold (DETECT_MINAREA) 
was set to be a minimum of 5 connected pixels with an intensity 
of0.5 5g (DETECT.THRESH) above the background lOesch et all 
|200'4 

lOesch et alj ( 120090 considered only those galaxies detected 
with S/N > 4.5 to be reliable. Since we were interested in galax- 
ies which are too faint to be individually studied, but that produce 
a relevant overall light contribution, we had to mask all the reliable 
single galaxies and focus on the total contribution of the fainter 
ones. We noted that, by doing this, we were sensitive to both galax- 
ies just below our detection threshold and much fainter ones. This 
was intentional, as we were trying to constrain all contributions be- 
low those of reliably detected galaxies. 

Starting from the segmentation maps, one of the outputs of 
SExtractor, we created two masks (one for each band) to reject all 
the detected sources existing in the field of view. Since most of 
the objects are surrounded by a bright halo, we decided to enlarge 
the masked area to avoid contaminating signal from bright galax- 
ies. We reversed each mask to be 1 where there are sources and 
elsewhere, we convolved two times the mask with a Gaussian filter 
with FWHM = 30 pixel = 2'.' 7 enlarging it each time by adding 
all the pixels with a value > 0.5. Finally we reversed again the 
mask. In this way we added ~ 18 — 19% more pixels to the mask. 

We obtained the final mask by merging the masks in the two 
passbands and convolving again. The masked area coiTesponds to 
16.8% of the field. 



3.2 The Power Spectrum Technique 

A description of random fluctuations is fundamental to extract in- 
formation hidden in the background of images and to provide a 
quantitative measure for the comparison with simulations. The spa- 
tial power spectrum is the key of our method, allowing us to deter- 
mine the amplitude of surface brightness fluctuations and to dis- 
tinguish the sources of noise from the fluctuations due to the unre- 
solved galaxie s. The power spe ctru m we mention is not the same 
applied first bv lHarrisonI ( Il970l) and lZeldovichI ( Il972l) . and then by 
[White ( 1 9940 to study the cosmic microwave background (CMB) 
anisotropies. We are, instead, interested in the spatial distribution of 
the faint z ^ 6 galaxies, so as to be able to extend the faint end of 
the luminosity function. For this reason we used the Fourier trans- 
form to derive the spatial power spectrum of the luminosity fluc tua- 
tion du e to these galaxie s , simi lar to the work of Tonrv & Schneideg 
( Il988l) and lTonrvetal.1 ( Il990l) . Using this method we are able to 



constrain the number of the faint galaxies (see, however, comments 
at the end of Section [33] l. 

A technique s i milar to ours was used for the first time by 
iTonrv & Schneideg (119881) . They discovered that the distance of a 
galaxy is inversely proportional to the amplitude of the luminos- 
ity fluctuations due to unresolved red giant stars and they used 
the spatial power spectrum to directly measure these fluctuations. 
The main difference between their work and ours is the density of 
sources per pixel. They were studying nearby galaxies that have 
a projected density of stars of many tens per pixel so the fluctu- 
ations between adjacent pixels have rms variations < 10% of the 
mean signal. On the contrary, we are dealing with images where the 
projected density of faint galaxies is much less then one per pixel, 
therefore we have smaller fluctuations and our results are more sen- 
sitive to the presence of spurious signals. 

Our approach consists in applying the IDLJj Fast Fourier 
Transform (FFT) routine to the masked images to compute the two- 
dimensional Fourier transform and, then, to derive the spatial power 
spectrum of the signal. Therefore, the power spectrum image ob- 
tained is two-dimensional, but to plot it and to do all the calcula- 
tions in the following it is convenient to derive the radial trend. 

The power spectra for the two bands can be compared by cal- 
culating their difference or ratio. The difference is the most natural 
choice, but relies on an extremely accurate color calibration of the 
power spectra in order to properly subtract the contribution of lower 
redshift objects from the zsso power spectrum. The ratio allows us 
to avoid the need of a very accurate and challenging color calibra- 
tion of the power spectra. Thus, from here on, we focus our analysis 
on the ratio of power spectra. 

We obtained the power spectra for both the J775 and Z850 im- 
ages (Figure [TJ. To highlight the light contribution coming from 
the undetected galaxies, which are bona fide 2: ~ 6 candidates, we 
examined the ratio between the zsso and 1775 power spectra (Fig- 
ure|2}. These undetected galaxies are responsible for a peak located 
between wavenumbers, k, of 100 and 400, shown in Figure|2] We 
tested the reliability of this feature using a x^ statistics. Comparing 
the values of the ratio between the power spectra in the range of 
the peak, i.e. 100 < fc < 400, and in the range 1400 < fc < 1700 
dominated only by the noise, we obtained xt ~ 3.08. This value, 
which can be rejected with a confidence well beyond 10a, con- 
vinced us that the peak is a real feature and it can not be ascribed 
to random noise. 

The angular scale 9, plotted on the top x-axis of the figures, is 
derived from k as follows: 



arcsec — 



Dframc [pixel] ■ scale [arcsec pixel ^] 
k 



(1) 



where Dframc is the dimension of the frame in pixels and scale 
is the pixel scale of the image. 

Following our technique, we can derive the number of faint 
galaxies from the height of the signal excess they produce in the 
ratio between the power spectra (Figure |2] and by comparing it 
with the one obtained from simulations (see Section|4j. 

It must be noted that our result can be affected by a frac- 
tion of interlopers at lower redshift. This is an issue affecting all 
studies based on the dropout selection and it is due to the aliasing 
betw een the Lyman break and the 4 000 A break, as discussed by 
Dahl en et alj <2010l) . ISu et al.l ( l201ir) . on the basis of accurate sim- 
ulations, estimated that the fraction of z ~ 1.2 galaxies which can 
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Figure 1. Power spectra, with colored en'or bars, of the background signal 
of the 4775 (top line with purple error bars) and zgso images (bottom line 
with blue error bars). 



Figure 3. Ratios between the zg^o and J775 power spectra obtained mask- 
ing all the sources detected by SExtractor with S/N > 4.5 (black line) and 
S/N > 3 (green dashed line). The percentages indicate the masked area. 
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Figure 2. Ratio between the zgso and 4775 power spectra. The shaded grey 
area represents the errors associated to the ratio. The region between fc = 
100 and k = 400, where the signal from the faint galaxies was measured, 
is highlighted in orange. 



be selected as 2 ~ 6 populatio n could be as high as 24%. This 
value is a pessimistic one since iMalhotra et al.l ( 120051) found half 
that number of contaminant galaxies, based on ACS grism obser- 
vations. 



3.3 Different Contributions 

The use of the power spectrum technique requires a preparatory 
study to understand how the different components of a scientific 
image affect the final power spectrum. 

The ratio between the spatial power spectra in Figure |2]can be 
divided into three different parts, the first one at low wave numbers 
(k < 100), the second one in the range of high wave numbers 
(k > 400) and the remaining one at intermediate wave numbers, 
which corresponds to the contribution from the undetected galaxies 
at 2 ~ 6. 

The power spectrum of random noise is a white power spec- 
trum having approximately a constant value over all the range of 
wave numbers. At high wave numbers, meaning at the pixel scale, 
the white noise contribution is dominant. At low wave numbers the 
principal contribution is due to residuals of the flat field correction 



process (for more details on this see Section 3.4). The highlighted 
region of the plot in Figure |2]is the most interesting one, where the 
signal is dominated by the contribution of the galaxies, therefore it 
permits us to infer the existence of galaxies below the usual detec- 
tion limit. In particular, the amplitude of the peak visible between 
k = 100 and k = 400 is proportional to the number of undetected 
galaxies in the field. 



3.4 Possible Sources of Contamination 

There are many sources of uncertainties that can affect our results 
and we will briefly list and characterize them. 

Mask: The primary source of uncertainty in our results is the 
size of the mask. Our aim is to mask all the galaxies which can be 
reliably and individually identified and to study the overall signal 
coming from those which are too faint to be detected one by one. 
If we are masking to a signal level which is too low, we will not 
see any signal coming from faint galaxies. On the other hand, if the 
mask brightness limit is too high the faint z ~ 6 undetected galax- 
ies will be overwhelmed by the light of the foreground galaxies and 
their bright halos. The ratio between the power spectra obtained by 
masking all the sources detected by SExtractor with S/N > 4.5 
and S/N > 3 are compared in Figure |3] In the region between 
k — 100 and k — 400 we notice that the peak due to the signal of 
faint i-dropouts is significantly lower and almost invisible. As the 
percentage of masked area is quite similar (16.8% and 20%), we 
can infer that we are detecting a signal from galaxies barely below 
our detection threshold rather than from a multitude of very faint 
objects. 

Calibration Files and Dither Pattern: The use of MultiDriz- 
zle could introduce spurious signals in the final scientific frame ob- 
tained after the combination of all the^f files. To test this, we down- 
loaded all the^f data available for the NICP12 field in the i775and 
085obands from the HST archive and replaced the scientific frame 
with random noise. The resulting mock frames have no sources and 
maintain the same coordinates of the real ones. Then, for each filter 
we processed these noise-only frames with MultiDrizzle, and we 
applied the power spectrum technique to the final scientific images. 
As it can be noticed from Figure |4] no features are introduced by 
MultiDrizzle neither in the single power spectra (top panel), nor in 
the ratio (bottom panel). Moreover, we did a similar test combining 
mock frames with noise multiplied by a random residual of flat field 
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Figure 4. Power spectra derived by combining tlie noise-only frames with 
MultiDrizzle (top panel) and their ratio (bottom panel). The flatness of the 
profile indicates that no features were introduced when combining the im- 
ages with MultiDrizzle. 
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Figure 5. Power spectra derived by combining with MultiDrizzle noise- 
only frames after being multiplied by a random residual of flat field with 
a maximum of 3% (top left panel) and 1% (top right panel) and their ratio 
(bottom panel). As in Figure|4] no spurious feature is introduced in the final 
ratio when combining the images with MultiDrizzle. 



(with a maximum of 0.3% and 1%) and, then, applying the power 
spectrum technique. It should be noted that the flat field residual 
in this case introduces a pixel-to-pixel uncertanty, not a large scale 
one. Once again, we concluded that MultiDrizzle is not creating 
spurious signals (see Figure|5} and, consequently, it is not affecting 
the final result of our technique. 

Cosmic Ray Residuals: The presence of cosmic ray residuals 
in the frame results in a gentle increasing slope in the range of 
high wave numbers, which depends on the intensity and number of 
pixels affected by the cosmic rays. 

Zodiacal Light: The zodi acal light is t he dominant source of 
sky brightness in the near-IR ( iReachlll997n . We were only inter- 
ested in the fluctuations due to the zodiacal light cloud chang- 
ing in the sky during the time of the observations. According to 
iKashlinskv et alj (2005), considering two observations of the same 
field taken six months apart and computing their difference, it is 
possible to isolate the zodiacal light fluctuation from the contri- 
bution coming from galactic a nd extragalactic objects. With this 
method IKashlinskv et alj <2005l ) were able to derive an upper limit 



for the zodiacal light contribution and to affirm that it did not af- 
fect their results on the cosmic infrared background fluctuations. 
Unfortunately, a few issues prevented us from applying the same 
analysis to the HST data. First of all, observations taken six months 
apart imply different orientations of the camera which could intro- 
duce detector-related systematics leaving a signature in the power 
spectrum. Moreover, ACS is affected by the contribution of scat- 
tered sunlight whenever the angle is close to (or less than) 90° and 
this could determine a wrong estimate of the fluctuation ascribable 
to zodiacal light. Finally, at the moment we do not have observa- 
tions of our fields taken with such a long temporal delay, so we 
cannot study relevant zodiacal fluctuations. On the basis of above, 
we decided to as s ume a negligible zodiacal contribution following 
IKashlinskv et alj 120051) and to better handle the zodiacal light fluc- 
tuation in our next paper, analyzing observations of the same field 
obtained with a long temporal separation, so as to highlight relevant 
differences in the zodiacal light. 

Flat Field Correction and Bias Subtraction Residuals: If the 
flat field correction or bias subtraction process were not correctly 
done, the data could be affected by a diffuse signal not coming 
from any real source. It can easily overcome the small amount of 
photons produced by the undetected galaxies we are looking for. 
These residuals modify the trend of each single power spectrum, 
increasing its slope in the low and medium wave numbers ranges 
and making difficult to see any peak in the power spectra ratio. 

Rotation: The images we used, oriented with North up 
and East left, were obtained from MULTIDRIZZLE without 
i ntroducing any additional correlation between adjacent pixels 
ICasertano et al. 2000) because of the use of the point kernel 
HBeckwith et al...2006. ). We did not introduce any further rotation 
in our data, nevertheless, we decided to investigate this effect for 
sake of completeness since this is a technical paper describing all 
the aspects of our new technique. If we were dealing with rotated 
frames it would be important to consider the following effect. The 
rotation of the images before extracting the power spectrum deter- 
mines a correlation between adjacent pixels, which is detected as a 
clear slope at high wave numbers of both the 1775 and zgso power 
spectra (Figure|6j. The magnified panel shows that the power spec- 
tra of the original and rotated data present different slopes, but, as 
shown in Figure|7] the ratio of the power spectra is not significantly 
affected by rotation since the slope introduced is exactly the same 
in both the 1775 and zgso power spectra. The ratio of the rotated 
data fits extremely well with that of the unrotated ones. 

Binning: To test the effect of data binning we compared 
the results obtained from the images with a pixel scale of 
90 mas pixel" ^ and a dimension of 3500 pixel each side, and the 
original images, characterized by a pixel scale of 30 mas pixel~^ 
and a 10500 pixel side. By binning the data, we lost details and so 
the power spectrum covers a more limited range of wave numbers. 
However, the ratio between the power spectra plotted in Figure |8] 
shows that the white power spectrum is well reproduced by both 
the 30 and the 90 mas pixeP^ cases. 

Groups of Galaxies: It is known that galaxies at high red- 
shift are not necessar ily randomly distributed, but can be grouped. 
lOverzier et al.l ( l2006h studied the clustering properties of galaxies 
at 2; ~ 6 finding an angular correlation for galaxies brighter than 
2850=28.5 mag on scales larger than 10" and no correlation on the 
same scales going down to 2:850=29. mag due to low S/N. Study- 
ing a sample of galaxies at z ~ 5, iLee et al.l ( I2OO6I) found a fur- 
ther correlation on smaller angular scales, which was interpreted 
as due to the so-called one-halo term. The clustering properties, 
due to both the one-halo and the two-halo term, are still not well 
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Figure 6. Power spectra obtained from the J775(top panel) and zggo (bottom 
panel) images with (bottom line with red error bars) and without (top line 
with blue error bars) applying a 45 ° clockwise rotation. In each panel the 
power spectrum referring to the rotation was vertically shifted to match the 
other one at fc = 125 (Dgg^i^jxies = 2'.'52) to better show the effect of the 
rotation. 



[arcsec] 



10 1 0.5 0.3 0.2 



0.1 



-0.8 




,0i,f00p0it0mimm i m mi^ 



v2.0 30 mas pix ' 
v2.0 90 mas pix"' 



1000 



2000 3000 
k 



4000 



5000 



Figure 8. Ratios between the zgso and J775 power spectra obtained for 
the v2.0 images with a pixel scale of 30 (blue solid line error bars) and 
90 mas pixcl"^ (red dashed line error bars). By biiming the data the power 
spectrum covers a smaller wave number range, but the different pixel scale 
is not affecting the features we are interested in. 
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Figure 7. Comparison of the ratios between the zssoand J775 power spectra 
for the v2.0 data with (red dashed eiTor bars) and without (blue solid error 
bars) applying a 45° clockwise rotation. The rotation has no significant 
effect on the ratio. 



explored for 1775 -dropouts with magnitudes in the range between 
2850=28 mag and zg5o=30 mag. Moreover, our field is too small 
to be affected by the large angular scale clustering which is bet- 
ter constrained JOverzier et al.ll2006l) than the small scale one. On 
the basis of this, we explored the effect of galaxy grouping on the 
power spectrum building up a toy model which introduces pairs of 
sources instead of randomly distributing them according to a non- 
zero angular correlation function only on scales smaller t han 10". 

W e assumed the one-halo term of uj{9) found by iLee et alj 
(120061) for Veoe -dropout galaxies to be valid also for 1775 -dropouts 
and placed the mock galaxies in the simulated frame accordingly. 
We found that clustering could amplify the signal from faint galax- 
ies. In particular, simulations with no clustering leads to an over- 
estimation of number of faint galaxies (see Section|4]for more de- 
tails). In any case, since the clustering/grouping properties of 2775- 
dropouts on the smaller angular scales are not known at the present 
epoch, our findings need to be revisited as soon as a proper angular 
correlation function is found at 2; ~ 6. 



4 MONTE CARLO SIMULATIONS 

The power spectrum permits us to see the existence of an overall 
signal coming from faint galaxies, but we are not able to directly 
obtain neither the number of those galaxies nor their magnitudes. 
To derive the number and properties of the undetected galaxies re- 
sponsible for the signal detected in the FFT we performed a series 
of Monte Carlo simulations. 

To create the mock 2775 and zsso images of the NICP12 field 
we started by deriving the values of the magnitudes, below the de- 
tection limit, and the morphology to assign to the simulated galax- 
ies. Then, we focused on the way to resemble as well as possible 
the noise and the large scale effects characteristic of the data. Fi- 
nally, we compared the power spectra obtained from the ACS data 
with the results coming from the simulations. We iterated the sim- 
ulation 100 times to estimate the variability of the power spectra 
depending on the elements we used. 



4.1 Elements of the Simulations 

In the following all the elements used in the simulations will be 
fully discussed. 

Random Noise: We used the IDL RANDOMN routine to cre- 
ate a frame with the same size of our images and randomly dis- 
tributed values with a mean of zero and a standard deviation of 
one. Subsequently, to reproduce the typical noise of the ACS data, 
we multiplied this frame by the characteristic standard deviation of 
the 1775 and 2850 images we measured with the IRAl|j task IM- 
STATISTICS. The resulting two random noise frames have a white 
power spectrum and reproduce exactly the trend observed at high 
wave numbers in Figure [T] 

Large Scale Effects: The residuals of the flat field correction 
produce the peak at low wave numbers in both the 1775 and 2350 



^ Imaging Reduction and Analysis Facilities (IRAF) is distributed by the 
National Optical Astronomy Observatories which are operated by the As- 
sociation of Universities for Research in Astronomy (AURA) under coop- 
erative agreement with the National Science Foundation. 
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Figure 9. Comparison of the ratios between the zgso and 1775 power spec- 
tra obtained from simulations including only noise and large scale residuals. 
Different models of large scale residuals were tested: a constant (green short 
dashed line), a linear one (magenta long dashed line), and a Gaussian one 
(orange dotted line). In black is plotted the ratio between the power spec- 
tra derived from the data. The dashed vertical line indicates the limit of the 
wave numbers range affected by the type of model considered for the large 
scale residuals. 



power spectra (Figure[T). Since we were not interested in studying 
tlie trend in tlie low wave number range, we simply searched a way 
to recover the shape of the power spectrum in this range. We created 
a model of flat field residuals to reproduce the large scale effects 
and, then, we added it to the simulated frame (additive flat field 
residuals). To be sure our model was not affecting the results of the 
simulations, we performed tests using mock frames with only noise 
and large scale residuals, assuming different models for the latter. 
In particular, we considered a Gaussian model, a uniform model 
(i.e. a constant flat), and one with a linear slope. As can be seen in 
Figure |9]the differences in the results are confined in the very low 
wave number range (k < 100), up to the dashed vertical line and it 
is not possible to reproduce the signal excess due to faint galaxies 
with only noise and flat field residuals. On the basis of these results 
we can affirm that the adopted model is not affecting the range we 
are interested in, and the signal in 100 < k < 400 is not affected 
by our choice regarding the modeling of the residuals of the flat 
field correction. 

Faint Galaxies: To determine the magnitudes of the simulated 
faint galaxies, we considered all the detected galaxies of the field. 
The histogram of the distribution of the values of the total magni- 
tude (obtained from the MAG_AUTO parameter in SExtractor) for 
the 2775 image is shown in the top panel of Figure [To] The value 
of 28.5 mag, where the number of galaxies decreases by a factor 
larger than 50%, is assumed as the brighter limit for the range of 
magnitudes used in the simulation of the J775 image. The fainter 
limit is 1775 = 30.5 mag, as we will see later in this Section. 

The number of faint galaxies with 2: ^ 5 to be inserted in the 
mock 1775 image was derived extrapolating the trend of the number 
counts beyond the cut off assuming a slope a; = — 1.3 according 
to [Metcalfe et al. ( 1995). The number of galaxies for each bin of 
magnitude used in the simulation is specified in Table |2] 

To reproduce the 2:850 image, we considered both the z ^ 6 
faint galaxies and those at 2: ^ 5. The latter include the same mock 
galaxies used to simulate the J775 image, whose magnitudes were 
corrected with a color term randomly drawn from a distribution re- 
producing the 1775-2850 color observed for the detected galaxies. 



600 

400 

200 


400 

200 




■ , -^-T^ 


r- 


r- 


r- 








-1 


1775 : 


■ ^^^^ 


p 


r- 












^850 



20 25 30 

total magnitude [ma^ 



35 



Figure 10. Distribution of the total magnitude of the detected galaxies in 
the J775 (top panel) and zgso (bottom panel) images. The vertical line cor- 
responds to the magnitude where we observe a drop in the galaxy number 
and which was assumed as the brightest magnitude used in our simulations. 



Table 2. Number and magnitudes of the simulated galaxies. 
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ai = -1.3 


2 = 6, 


a= -1.9 
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Number 


28.5 


714 


28.0 


32 


29.0 
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28.5 


87 


29.5 
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29.0 
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30.0 


5503 


29.5 


622 


30.5 


10028 


30.0 


1665 



The former population consists of 2: ~ 6 galaxies not detectable 
in the 1775 image. The brighter limit for the magnitude range is 
assumed to be the value of 28 mag where the histogram of the to- 
tal magnitudes shows a drop larger than 50% (Figure [lO] bottom 
panel). The fainter limit of the range considered is 2850 = 30.5 
mag. We extrapolated to fainter magnitudes the luminosity func- 
tion at z ~ 6 to obtain the number of galaxies to be created. 

All the simulations hereafter include a flat field residual 
model. After a series of tests, we found a good agreement in re- 
producing this effect by combining 20 Gaussian functions with a 
standard deviation of about one third of the image dimension: 10 
of them have a peak of 10^* to match the typical fluctuations of 
the flat field, and the other 10 Gaussians have a peak of —10^* to 
create a smooth background with soft large scale fluctuations. The 
chi-square statistics results obtained comparing the data to the sim- 
ulations including the three different flat-field models plus mock 
galaxies and noise in the range 100 < k < 400 were very similar 
(see Columns 3 and 4 in Table[3]and panel (b) of Figure II 6t. 

We tested different values for the slope a, in particular a — 
-1.3, a = -1.7, Q = -1.9, a = -2.1, and a = -2.4. On 
the basis of a x^ statistics (see columns 7 and 8 of Table |3] and 
(b) panel of Figure 116) we assumed our extrapolation to have the 
slope a — —1.9 (Figure [TTJ and we anchored it to t he LF value 
at a magnitude of 28 mag given bv iBouwens et al. ( 12007 ). This 
value is in agreement with the expectation of lSu et al.l ( 1201 ih . who 
found the faint end slope of the LF at z ~ 6 to be in the range 
-1.90 < a < -1.55. 

Testing the effect of groups of galaxies as described in Section 
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Figure 11. Comparison of the ratio between the simulated power spectra 
derived assuming different faint end slopes. The green, red, orange, blue, 
and magenta lines represent the median obtained from simulations assum- 
ing a = —1.3, a = —1.7, a = —1.9, a = —2.1, and a = —2.4, 
respectively. In black the ratio derived from the data is plotted with its error 
bars. 
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Figure 12. Comparison of the ratio between the simulated power spectra 
derived assuming different faint end slopes and pairs of galaxi es with a 
distan ce derived from the angular correlation function at ^ ~ 5 iLee et alj 
|200g). The green, red, and blue lines represent the median obtained from 
simulations assuming a = —1.7, a = —1.8, and a = —1.9, respectively. 
In black the ratio derived from the data is plotted with its error bars. 



13.41 we considered different slopes a, i.e. a — —1.7, a = —1.8, 
and a — —1.9 introducing pairs of mock galaxies in the simulation 
(Figure[T2j- The x^ statistics shows that a = —1.8 fits best the data 
(xl = 1.06, see Columns 9 and 10 of TablelDl. 

Finally we also considered the effect of contaminants at lower 
redshift changing the low-luminosity slope of the distribution of 
galaxies with 2 ^ 5. To this aim we considered higher and lower 
values for Oi and we noted that this change produces a similar one, 
with opposite sign, in the derived a value at 2; '^ 6. 

In building the set of galaxies, we considered a range of ge- 
ometrical (i.e., position and orientation on the frame) and struc- 
tural parameters (i.e., effective radius, r^, and axis ratio, b/a). Each 
synthetic galaxy was generated with a set of random, equally dis- 
tributed, values for the former parameters and a set of values which 
resembles the characteristics of the fainter detected galaxies of the 
NICP12 field, derived from the catalog obtained running SExtrac- 
tor, for the latter. Regarding the former parameters, each synthetic 
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Figure 13. Comparison between the ratio of the power spectra obtained 
from the data (black line) and from the simulations assuming our effective 
radius distribution (magenta line) and the one by Ferguson et al. (2004, 
green line) 



galaxy was generated with a set of random values. On the other 
hand, the structural parameters values had to resemble the charac- 
teristics of the fainter detected galaxies of the NICP12 field, de- 
rived running SExtractor. The values of the effective radius were 
chosen within the range 1.5 ^ re $J 11 pixel corresponding to 
0'.'135 ^ ro ^ Cl.'99, and the range spanned by the axis ratio 
was 0.3 ^ b/a ^ 1. Even though the distribution of re we used 
was o nly slightly different from the one found by IFerguson et al.l 
( 12004) at z ~ 5, we tested the effect of this parameter on the sim- 
ulations (Figure [Tit. In particular, the reduced chi-square statistics 



(xt 



1.22) supports the idea that small differences in the Tc dis- 



tribution do not significantly affect the result. 

Finally, we convolved the moc k galaxies with the HUDF PSF 
(FWHM ~ rt'l, IOeschetal.l2007h . 

Galaxy Morphology: Regarding the morphology of the mock 
galaxies, we adopted a Sersic profile and, changing the Sersic in- 
dex n, we obtained ellipticals (n = 4) and spirals (n — 1). Since 
there are no clear indications regarding the distribution of different 
galactic types in the early Universe, we studied the effect of chang- 
ing the relative number of ellipticals and spirals in our simulations. 
We compared the results obtained introducing 15% ellipticals (de- 
scribed by a de Vaucouleurs profile) and consequently 85% spirals 
(described by an exponential profile) in the simulations with those 
obtained simulating 50% ellipticals and 50% spirals and 33% ellip- 
ticals and 67% spirals. In Figure [14] the ratio of the power spectra 
obtained from the data is compared with the mean ratios and 3cr 
confidence ranges derived from simulations with different distribu- 
tions of ellipticals and spirals. Applying the chi-square statistics to 
the wave number range where the peak is, i.e. 100 < fc < 400 (see 
Columns 1 and 2 in Table[3} we selected the realization character- 
ized by the smallest x^ value as the best fit (see also the distribu- 
tion of the residuals in the (a) panel of Figure 1 16). It corresponds to 
the mix of 33% ellipticals and 67% spirals. These perc entages are 
consistent with the findings of ISomerville et al.l ( 12001 ), who pre- 
di cted the majority of ga laxies at 2 ^ 3 to be disk dominated, and 
of iLotz et all ( I2OO6I) and lRavindranath et all ( l2006h who found the 
percentage of galaxies with bulgelike morphology to be ~ 30% at 
2 ~ 3 and 2 ~ 4, respectively. 

This mix was adopted to test the sensitivity of the method to 
the assumed value of the faint end limit. To this aim we performed 
a few simulations changing the fainter magnitude limit of the mock 
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Table 3. The results obtained from the x^ statistics are listed for all the realizations of the simulations regarding the percentage of elliptical and spiral galaxies 
(Columns 1 and 2), flat field residual model (Columns 3 and 4) introduced, value of the cut-off in magnitude (Columns 5 and 6), and slope of the faint end of 
the LF with single sources (Columns 7 and 8) and with pairs of galaxies (Columns 9 and 10), respectively. The selected realization is the one characterized by 
the smallest xi value, i.e. xi = 1-02. A value of xt = 1-24 con'esponds to a probability of 0.3% and it can be rejected with a 3(t confidence level. 
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Figure 14. Comparison of the ratio between the zgso and 1775 power spec- 
tra (black line) with the median derived from the simulations with a different 
percentage of spirals and ellipticals represented by different colors, 50% el- 
lipticals and 50% spirals in green, 33% ellipticals and 67% spirals in red, 
and 15% ellipticals and 85% spirals in blue. 



[arcsec] 



1 0.7 




100 



200 



300 



400 



Figure 15. Comparison of the observed ratio between the zgso and 1775 
power spectra (black) and those obtained by simulating unresolved galaxies 
in different magnitude ranges, i.e., 28.5 < 1775 < 29.5 (blue dashed fine), 
28.5 < 1775 < 30 (orange dotted line), 28.5 < 1775 < 30.5 (red dashed 
line), and 28.5 < 1775 < 31 (green dotted line). 



galaxies (Figure ITSt. We used the chi-square statistics to select the 
best simulation concerning the magnitude cut-off (see Columns 5 
and 6 in Table |3). Moreover, we computed the difference between 
the simulations and the data, and we plotted the residuals (Figure 
IT6l (c) panel ). On the basis of above, we decided to create galaxies 
with magnitudes up to 30.5 mag. 

The Mask: Consistent results can be achieved only by ana- 
lyzing the same area in data and simulations. For this reason we 
used, for the simulated images, the same mask created from the 
ACS data. Moreover, we considered what percentage of the frame 
is masked and calculated the number of mock galaxies that had to 
be outside the masked area. We, therefore, placed them randomly 
on the frame checking that they were positioned outside the mask. 

The sum of the simulated galaxies, the frame reproducing the 
large scale effects, and the random noise frames was multiplied by 
the mask and analyzed with the power spectrum technique. 



4.2 Results of the Simulations 

The results obtained by iterating the simulations 100 times consid- 
ering the 5 bins of magnitudes given in Table |2] and assuming one 
third of the galaxies to be ellipticals and the others spirals are shown 
in FigurefTT] The trend of the data with the associated errors is well 
fitted by the mean ratio between the simulated power spectra, and 
all the values are in agreement with the 3a confidence range. 



4.3 The z ~ 6 Faint Galaxies Contribution to Cosmic 
Reionization 

In order to evaluate the implications for reionization of our findings, 
we considered the total surface brightness contribution of the un - 
detected galaxies following the approach of lStiavelli et aU ( l2004ah . 
When this contribution is added to that of detected galaxies, we find 
a surface brightness of 24.8 mag arcmin^^ when extrapolating to 
a magnitude zsso = 29 mag and 23.6 mag arcmin"^ when extrap- 
olating to 2850 = 30 mag. We compared these values to the min- 
imum surface brightness of 27. 2 mag arcmin~^ need ed for reion- 
ization by Population II stars Jstiavelli et alj|2004bf) . Our results 
give comfortable margins of a factor 9 and 27, respectively, to ac- 
commodate the escape fraction /c and the clumping factor C. For 
instance, assuming C = 6 and the Populat ion II spectral energy 
distribution (SED) of IStiavelli et all ( l2004bh . the LF extrapolated 
down to 2850 = 30 produces enough ionizing photons to reionize 
the Universe if /c > 5%. 

So far the implications for reionization of what we found seem 
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Figure 16. Distribution of the residuals derived computing the difference 
between the different realization of the simulation and the data, in terms of 
the ratio between power spectra in the range between k=100 and k=400. 
The top left panel, (a), refers to the different percentages of spirals and el- 
lipticals, the centi'al left one, (b), to the various flat field residual models, the 
bottom left one, (c), concerns the magnitude cut-off, and the right one, (d), 
is about the different value of the faint end slope a. The filled histograms 
refer to the selected realization, i. e. the one with the lowest x^ ■ 
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Figure 17. Mean ratio between the J775 and zgso simulated power spectra 
(blue line) and Scr confidence range (light blue region) derived from 100 
simulations compared to the results coming from the ACS data (black hne). 



foitn to obtain the spatial power spectram of the background signal 
and we examined the ratio between the zsso and 1775 power spectra. 
In this way we were able to distinguish the sources of noise from 
the fluctuations due to the unresolved galaxies. We tested the effect 
of several parameters on the power spectrum to test its reliability 
and how contamination can affect the results. 

We showed that the power spectrum technique is a powerful 
tool in analyzing the light contribution produced by galaxies which 
are below the detection limit in deep and ultra deep surveys. We 
used this technique to estimate the contribution to cosmic reion- 
ization from faint galaxies (zsso ^ 28 mag), which are bona fide 
2 ~ 6 candidates, in the NICP12 ACS field of the UDF05. 

Monte Carlo simulations were used to determine the num- 
ber of faint undetected i-dropouts responsible for the peak ob- 
served in the ratio between the zsso and 1775 power spectra (Fig- 
ure |2}. The data are consistent with a faint end slope of the lu- 
minosity function of a = —1.9 if not introducing any cluster- 
ing/grouping in the 1775 -dropouts distribution, and a — —1.8 
when considering pairs of galaxies with distances drawn accord- 
ingly to the angular correlation function by Lee et al. (2006). Both 
th e Q values we ob t ained are consistent with p revious findings 
bvlYan & Windhorstl(l2004t).lMalhotra et al.lll2005l),lBouwens et al 



2006l),lBouwens et al.1 (l2007h.lBouwens et alj J201 ih. andlSu et al 



2011D. 



Considering a = —1.9, adding the total surface brightness 
contribution of faint galaxies to that from the detected galaxies and 
comparing it to the minimum value required to reionize the Uni- 
verse, we obtained a margin to model the escape fraction and the 
clumping fact or. Adopting the clum ping factor of 6 and Popula- 
tion II SED of lStiavelli et all ( l2004bl) . the z ~ 6 undetected galax- 
ies down to 2850 = 30 mag produce enough ionizing photons to 
reionize the Universe assuming the escape fraction to be larger than 
~ 5% The solution with all galaxies in pairs and a slope a = — 1.8 
changes the cumulative surface brightness by less than 0. 1 mag and 
implies an escape fraction larger by ^ 10%. Thus, the main conclu- 
sion of the paper that reionization is indeed consistent with being 
completed due to the contribution of faint galaxies is not affected 
by the presence of multiple halo occupation. 

We expect to apply our technique to the ultra deep data ob- 
tained with NICMOS and the Wide Field Camera 3 (WFC3). The 
comparison of the results for different fields and instruments will 
permit us to better constrain the role that faint galaxies at 2: ~ 6, 
2 ~ 7, and z ~ 8 had in the cosmic reionization. 



to be very promising, but we are waiting for a more thorough anal- 
ysis of all available fields before pushing this interpretation any 
further. 



5 CONCLUSIONS 

In this paper we presented a new technique to quantify the light 
contribution coming from the faint high redshift galaxies below the 
detection threshold of imaging data. Specifically, we illustrated the 
technique with an application to ACS data obtained in the F775W 
and F850LP filters for the HUDF parallel field NICPI2 looking 
for faint 2: ~ 6 galaxies on the basis of the Lyman break technique, 
which is based on a comparison between the flux of a galaxy within 
the set of two filters. 

In particular, our aim was to extend by a few magnitudes the 
faint end of the luminosity function at that redshift. After masking 
all the detected sources in the field we applied a Fast Fourier Trans- 
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